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Abstract 

We analyze the Lagrangian flow in a family of simple Gaussian scale- 
invariant velocity ensembles that exhibit both spatial roughness and tempo- 
ral correlations. We show that the behavior of the Lagrangian dispersion of 
pairs of fluid particles in such models is determined by the scale dependence 
of the ratio between the correlation time of velocity differences and the eddy 
turnover time. For a non-trivial scale dependence, the asymptotic regimes 
of the dispersion at small and large scales are described by the models with 
either rapidly decorrelating or frozen velocities. In contrast to the decorre- 
lated case, known as the Kraichnan model and exhibiting Lagrangian flows 
with deterministic or stochastic trajectories, fast separating or trapped to- 
gether, the frozen model is poorly understood. We examine the pair disper- 
sion behavior in its simplest, one-dimensional version, reinforcing analytic 
arguments by numerical analysis. The collected information about the pair 
dispersion statistics in the limiting models allows to partially predict the 
extent of different phases and the scaling properties of the Lagrangian flow 
in the model with time-correlated velocities. 



1. INTRODUCTION 

The aim of this paper is to study the Lagrangian flow in (i-dimensional random velocity 
fields v(t,r) with a prescribed scale-invariant statistics. The velocity ensembles that we 
shall consider mimic some essential properties of realistic velocities in developed turbu- 
lence: their spatial roughness within a large interval of scales and their temporal correla- 
tion. By definition, the Lagrangian flow is described by the ordinary differential equation 
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(1.1) 



It determines the motion of hypothetical fluid particles or of small test particles suspended 
in the fluid. One usually distinguishes between the motion of a single particle, dominated 
by the velocity fluctuations on the largest scale present (the so called "sweeping effects") 
and the evolution of a relative separation of two particles, the object of interest of the 
present paper, driven by the velocity fluctuations on scales of the order of the inter-particle 
distance. For larger groups of particles, one should similarly distinguish the motion of their 
barycenter from the relative motion of particles within the group. The latter is known to 
show quite intricate behavior related to intermittency, see [10], but it will not interest us 
here. 

The separation p = R — R between two fluid particles satisfies the equation 

^ = v(t,p + R(t))-v(t,R(t)), (1.2) 

where R(t) is a trajectory of one of the particles, a solution of Eq. (1.1) starting at time 
zero at R = 0, for example. Upon introduction of the so called quasi-Lagrangian velocity, 

v qL (t,r) =v(t,r + R(t)), (1.3) 

i.e. velocity in the frame moving with a fixed fluid particle, we may rewrite Eq. (1.2) as 

^ = v« L (t, p) - v« L (t, 0) = Av qL (t, p) . (1.4) 

We shall be interested in the short- and long-time behaviors of the Lagrangian particles in 
the statistical ensembles where typical velocities are only Holder continuous in space. In 
such non-Lipschitz velocities, there is a problem with solving Eqs. (1.1), (1-2) or (1.4). To 
avoid it, we shall first consider noisy particle trajectories that solve the stochastic equation 

dR = v{t,R)dt + \PlkdW , (1.5) 

where W(t) is the Brownian motion in d dimensions. Noisy trajectories form a well 
defined Markov process even in velocity fields with poor regularity. Subsequently, the 
limit k, — > will be performed in selected quantities. The velocity ensembles that we 
shall consider are time- reversal invariant. As a result, we shall not have to distinguish the 
forward and the backward evolution of trajectories and will concentrate on the first one, 
between, say, times and t. 

One way to study the relative motion of pairs of fluid particles is to follow the evolution 
of the pair dispersion, i.e. of the separation distance p between two particles. Its statis- 
tics in a random flow may be described by the velocity-averaged probability distribution 
V(po, p',t) dp of the time t dispersion p given its time zero value p , in the limit when 
we remove the (independent) noises of the Lagrangian trajectories. As we shall see, the 
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PDF's V(po,p;t) may be sometimes distributions rather than integrable functions. An- 
other, related, test of the relative motion of a pair of fluid particles is obtained by looking 
at the exit time [17,4]: the time t that the pair dispersion takes to evolve from p to p x . 
In particular, p\ = 2/?o corresponds to the doubling time of the pair dispersion. The 
statistics of the exit times may be encoded in their velocity- averaged probability distribu- 
tion Q(t] Po, Pi) dt taken in the limit of vanishing noise. The exit time is less influenced 
than the pair dispersion by small or large-distance cutoffs in the velocity correlations, so 
preferable in numerical or experimental studies [4]. 

Recently, a new insight into the intricate character of the Lagrangian flow in turbulent 
velocities has been gained, see [3,19,24,30,25], by analytic study of the Kraichnan ensemble 
[23] of Gaussian velocities which are decorrelated in time but exhibit scaling behavior in 
space. Here, we try to find out how the presence of temporal correlations of velocities 
influences the Lagrangian flow. We shall study the behavior of trajectory separation in a 
simple generalization of the Kraichnan ensemble of velocities where time correlations are 
reintroduced. 

The paper is organized as follows. In Sect. 2 we recall the main facts about the La- 
grangian flow in the Kraichnan model, in particular the appearance of phases with very 
different trajectory behavior. Sect. 3 describes a Gaussian ensemble of velocities with tem- 
poral correlations, discussed in the past in [8,2,11,14]. It exposes a simple mean-field type 
analysis of the particle separation when such an ensemble is used to model the quasi- 
Lagrangian velocities. How the mean-field predictions may be substantiated further by 
scaling arguments is the subject of Sect. 4. See also [12,13] for related rigorous results. An- 
alytic arguments and conjectures about the behavior of trajectories in a one-dimensional 
version of the model with time- independent velocities are contained in Sect. 5. The behav- 
ior of the exit time statistics in velocity ensembles with long-time correlations is studied 
in Sect. 6. The question how the behavior of pair dispersion changes when the Gaussian 
ensemble is used to model the Eulerian velocities is addressed in Sect. 7. In particular we 
show that in the one-dimensional time-independent case, the sweeping by large eddies in 
the Eulerian model speeds up the movement of a single Lagrangian particle but it localizes 
pairs of particles by reducing the growth of their dispersion. Three Appendices contain 
more technical material. Some of the predictions of the paper are checked in one dimension 
by numerical simulations. 
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e a Tecnologia grant: PRAXIS XXI/BD/21413/99. KG, PH, AK and MV are grateful 
to the Erwin Schrodinger Institute in Vienna and the Institute for Advanced Study in 
Princeton where parts of the work on the paper were done. KG acknowledges the support 
of the Neumann Fund and MV the Ralph E. and Doris M. Hansmann Membership at 
the IAS. AK and MV were also supported by the European Union under the Contracts 
FMRX-CT98-0175 and HPRN-CT-2000-00162, respectively. 
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2. LESSONS FROM THE KRAICHNAN MODEL 



The Kraichnan ensemble of turbulent velocities [23], is a Gaussian ensemble with van- 
ishing velocity 1-point function and with the 2-point function 

/ \ rp ik-(r-r') Jfc 

(« f (*, rMf.rO) = D, 5(t-t>) J P-(k,p) — , (2.1) 

where k L = Vk 2 + L~ 2 and P ij (k, p) = ^(S ij - + p^f, see [10] and references 
therein. There are two dimensionless parameters in the Kraichnan ensemble: the exponent 
£ > and the compressibility degree < p < 1. For £ < 2, the velocities smeared in time 
are (almost surely) Holder continuous in space with any exponent smaller than £/2. For 
£ > 2, they are Lipschitz (or even more regular). The compressibility degree p = cor- 
responds to incompressible velocities, p — 1 to gradients of a potential (in one dimension 
necessarily p — 1). The normalization constant -Di has dimension length 2 ~^ /time. The 
length L is the integral scale that sets the spatial correlation length of velocities. If £ < 2, 
it may be taken to infinity in the correlation functions involving only velocity differences 
v(t, r + p) — v(t, r) = Av(t, p). For £ > 2, this may still be done if D x is rescaled when 
L — > oo. 

A. Possible flow behaviors 

The Kraichnan ensemble may be used invariably to model Eulerian or quasi-Lagrangian 
velocities as both ensembles coincide in this case. The statistics of a single Lagrangian 
particle is that of a rf-dimensional Brownian motion with diffusivity that blows up when 
L — > 0. The two-particle separation depends only on velocity differences and its statistics 
has a non-trivial L — > oo limit whenever this holds for the velocity differences. The pair- 
dispersion and the exit time PDF's V(po,p;t) and Q(t; po, p±) may be solved analytically 
in this limit. The exact solutions, that describes also the short-distance asymptotics of 
the PDF's at finite L, show several dichotomic behaviors depending on the values of 
parameters of the model. The first dichotomy, noticed in [3], is between the 

deterministic flow characterized by the property 

lim V(po,p;t) = 5(p) (2.2) 

po-*0 

which signals that the trajectories in a fixed velocity field are defined by their initial 
position, and the 

stochastic flow where 

lim V(po,p;t) is an integrable function of p. (2.3) 
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The limits of the PDF's above (and below) should be understood in weak sense, under 
integrals against test functions. The behavior (2.3) means that infinitesimally close tra- 
jectories separate in a finite time and indicates that the stochasticity introduced into the 
Lagrangian flow by coupling it to the noise, see Eq. (1.5), survives in the limit k — > 0. 
The Lagrangian trajectories in a fixed velocity field are not determined by initial position 
but form instead a stochastic process. That this is indeed what happens in the Kraichnan 
model was established rigorously in [24]. We shall call the phenomenon spontaneous 
stochasticity 1 . 

There are further dichotomic behaviors of the Lagrangian flow in the Kraichnan model. 
We have chosen to characterize the other dichotomies in terms of the small po behavior 
of the exit time PDF Q(t; p , 7p ) with < 7 = const. 7^ 1, attaching to them more or 
less suggestive names. There is a dichotomy between the 

Lyapunov flow such that 

lim <2(t;po,7Po) is an integrable function of t, (2.4) 

po— »o 

and the 

Richardson flow where 

lim Q(t;p ,7Po) = c( 7 ) S(t) with c( 7 ) > . (2.5) 

po— >0 

This dichotomy distinguishes the flows in regular velocities where the Lagrangian sepa- 
ration on short distances involves fixed time scales (like the inverse Lyapunov exponent), 
from the ones in non-regular (non-Lipschitz) velocities where the characteristic times of 
the Lagrangian separation become very short on short scales. 

Finally, the last two dichotomies that we want to single out characterize the short 
distance behavior of the probability of infinite exit times. They are between the 

locally separating and locally trapping flow where for 7 > 1 

Jim J Q(t;p ,-fp )dt { < } 1 , respectively, (2.6) 

and 

locally recurrent and locally transient flow where the same holds for 7 < 1. 

Roughly, with positive probability, close trajectories do not increase their distance in 
locally trapping flows and do not approach each other in locally transient ones. 



It was termed intrinsic stochasticity in [30] 
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B. Pair dispersion 

Due to the temporal decorrelation of the Kraichnan velocities, the probability measures 
V(po, p; t) dp constitute transition probabilities of a Markov process p(t) that in the limit 
L — > oo becomes a diffusion on a half-line with the explicit generator 

M = -D[ pt- a d p p a d p for a = £± - 1 , (2.7) 

where D[ is proportional to D\. Note that the symbol of M vanishes at p = 0. The 
different behavior of such diffusion for different values of £ and a has its origin in the 
singularity of M at p = which requires different treatment in different regimes. Up to 
the change of variables x = p^^l 2 casting the generator M into the form 

M = \(2-i) 2 D>x 1 - s d x x s - x d x for 5 = 2(1 - ±=|) (2.8) 

and a time rescaling, the Markov process p(t) may be identified with the well studied 
Bessel diffusion [5] , a natural interpolation between processes describing the radial variable 
in the standard diffusion in dimensions 5. Various analytic formulae may be then directly 
carried from that case to the present situation. 

For £ = 2 corresponding to smooth velocities with velocity differences linear in space, 
the particle dispersion PDF takes a log-normal form [7,9]: 

V(p ,p;t) oc e 4I V p , (2.9) 

where A = y^f^-D'i is the (biggest) Lyapunov exponent. It is easy to see that in this 
case, lim V(po,p;t) = 5(p), pointing to the deterministic nature of the flow. Indeed, in 

po^O 

velocities regular in space the trajectories are uniquely determined by their initial position 
and very close fluid particles separate little in a fixed time interval. Nevertheless, all 
moments of the pair dispersion behave exponentially in time and grow in the chaotic 
regime where A > 0, i.e. p < d/4, whereas sufficiently small (fractional) ones decrease 
when A < 0, i.e. p > d/4. For the second moment, one obtains: 

(p 2 (t)) = e ( 2A+4D ^Po- (2-10) 

Similar behaviors persist at small times in the Kraichnan velocities with £ > 2 and finite 
L, see [24]. 

The stochastic Lagrangian flow occurs in the non-Lipschitz version with < £ < 2 
of the Kraichnan model for weak compressibility p < d/C, 2 . In terms of the parameter 
b = = that will be frequently used below, the latter inequality means that b < 1. 
In this region, V(po,P',t) = e~ tM (po, p) where M is taken with "singular Neumann" or 
reflecting boundary condition at p — 0, see [19], and 
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lim V(0,p;t) oc p^t 1 " 1 e ^ D 'S . (2.11) 

Note the stretched-exponential form of the PDF. In particular, one obtains for the second 
moment of the separation distance and large times: 

(p 2 (t)) = 0{t^i). (2.12) 

This is the Kraichnan model version of the 1926 Richardson law [28] stating that in 
developed turbulence the mean square dispersion grows like t 3 . Note that the Richardson 
behavior is reproduced in the Kraichnan model for £ = |. In the limit when the initial 
trajectory separation p — * 0, the power law behavior (2.12) extends to the entire time 

domain and (p 2 (t)) oc t~ . 

In the non-Lipschitz strongly compressible version of the Kraichnan model correspond- 
ing to < i < 2 and p > d/£ 2 (or 6 > 1), 

V( Po ,p;t) = V^(p ,p;t) + p( Po ;t)S(p) (2.13) 

with the regular part of the PDF V reg (po, P',t) = e~ tM (po, p) where M is taken with 
"singular Dirichlet" or absorbing boundary condition at p — 0, see [19], and with 

p(po; t) = l - 7 (b, (^[i) r(rl ' (2 - 14) 

X 

where 7(6, x) = / e~ y ' dy is the incomplete gamma-function. When p — > then the 



regular part of the PDF tends to 0, whereas p(po; t) tends to 1. We recover this way the 
deterministic behavior (2.2), with trajectories in fixed velocity realizations determined by 
their initial position. The presence of the term proportional to S(p) at finite po signals 
that the trajectories starting at different initial positions coalesce with positive probability 
[19]. The time growth of the mean distance square dispersion is different here: 

(p\t)) = 0(pl-U^- b ) (2.15) 
with a logarithmic correction at p — <i/£ 2 i.e. at 6=1. 

C. Exit time 

In the Kraichnan model, the exit time PDF Q(t; po, pi) may also be directly controlled 
using the kernels e~ tM °(po, p) where denotes the generator M of (2.7) with the 
Dirichlet boundary condition at pi (in addition to the appropriate condition at the origin 
when pi > po). This is due to the Markov property of the stochastic process p(t). Since 
the exit times have not been discussed in the context of the Kraichanan model, we shall 
use the occasion to provide more details, essentially translated form the Bessel diffusion 
case. The PDF of the time of exit through pi is given by 
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Q(t; P o,Pi) = -d n(pi) e- tM °(p ,p), (2.16) 

where d n (pi) = ± D[ p a d p p^~ a \ p=Pl plays here the same role as the normal derivative in 
the classical potential theory. The sign is that of (pi — po). This PDF does not have to 
integrate to 1. The eventual missing probability corresponds to events where p(t) stays 
forever in the open interval (0,pi) or (p 1 , oo) or when it gets absorbed at the origin, see 
below. We shall assign the infinite value of the exit time to such events. The averages of 
the powers of the exit time t over the realizations with t < oo may be expressed by the 
kernels of the inverse powers of operator M D : 

oo 

(t n l { t<oo } ) = Jt n Q(t; Po , Pl )dt = -nld^M^-^p-t), (2.17) 
o 

where by 1{a} we denote the characteristic functions of the events satisfying the condition 
A. In particular, the probability that the exit time is finite 

oo 

(1{*<oo}> = / Q(t;p ,Pi)dt = -d n[pi) M D \p ,p;t). (2.18) 
o 

The expectations (2.17) may be obtained from the characteristic function 

<e^l {t <oo}> = -d n(pi) (M D -iu;)- 1 (p ,p 1 ;t) (2.19) 

involving the resolvent kernel of Mp. We may also consider the averages conditioned on 
the exit times being finite: 

\L{t<oo}) \l{t<oo}) 



For £ = 2 and L = oo, i.e. in the smooth version of the Kraichnan model, the kernel 
e~ tM ° (p , p) is easily calculable by the image method: 

e- tMD ( Po ,p) = 1 ( C -^ (HP/P0) - Xt)2 c ^ ln (pi/po)- i ^ 7 (ln(ppo/pD-A f ) 2 s ^ 

For the exit time PDF, one obtains: 

= \Hpi/po)\ c ^f a2 c -i^ 1 ° 2 ^/po)-^ (2 22) 

Note that the PDF depends only on ^ so that the corresponding Lagrangian flow is 
Lyapunov in our terminology, see (2.4). The total mass is given by the simple expression: 
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o 



if A ln( Pl /p ) > 0, 
Q{1:Ptt -" l)d1 < — if Aln( Pl /Po)<0. (2 ' 23) 




The missing probability corresponds for p\ > p and the negative Lyapunov exponent to 
the events when p(t) stays forever in the open interval (0, pi) (the pairs of trajectories 
remain always close). For p\ < p and the positive Lyapunov exponent, it represents 
the events when p(t) stays forever in the interval (p 1; oo) (the pairs of trajectories never 
come close). According to the characterization from the previous subsection, see conditions 
(2.6), the Lagrangian flow is locally separating and locally transient if A > 0, i.e. p < d/A 
and it is locally trapping and locally recurrent if A < 0, i.e. p > d/A. Finally, when 
A = 0, i.e. p — d/A, it is locally separating and locally recurrent. Similar properties of 
the exit time PDF hold for £ > 2 and finite L asymptotically at short distances. 

The exit time PDF (2.22) vanishes with all derivatives at t — and it has an expo- 
nentially decaying tail at large t. For the exit time conditional moments, one obtains: 



l,n\ _ |A| / |ln(pi/pp)| \ n+ ? 2D' /P0) ' ts- ( I A ln(pi/pp)| \ 



(|„_i|_i +fe )! f D[Y f\HPi/Po)\ 



fe! (| n _i|_i_ fc)! yx?) \ 



(2.24) 



where the first expression with the Bessel function holds for all real n and the second one 
for integer n. The unconditioned moments diverge for n > if A ln(pi/p ) < due to 
the finite probability of infinite exit times. For the conditional characteristic function one 
obtains: 



where the square root is taken with the positive real part and the sign is that of A ln(pi/p )- 
The decay at large |a>| and the analyticity properties of the characteristic function reflect 
the behavior of the exit time PDF at small and large t. In particular, the decay ~ 
e ~l Mpi/po)\\ / \^\/d 1 q £ charac^erig^ic function along the positive imaginary axis of uj 

signals the behavior ~ e 4D i* ( - pl ^ p0 " > f Q(t; p , pi) for small t and the singularity at 



uj = —i^y indicates the presence of the tail ~ e 4D i for large t, in agreement with 
Eq. (2.22). Note non-Gaussian character of the large fluctuations of t and of t -1 . 

In the non-Lipschitz version of the Kraichnan model with < £ < 2, the resolvent 
kernel (Mo — iu)^ 1 (p , p) may still be easily calculated in a closed form. Let us start 
from the case with p\ > p . On the interval [0,pi], 

(M D -iu)-\p G ,p) = Iw-V^/tU)" 1 



x 2 

4D 7 : 



-t 
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Mpo) U±(p) Mpi) - Mp) /±(Pi)) for Po < p, (9 9(]] 

U(p) (f±(po) U(pi) - U(po)f±(pi)) for Po >P, 1 ' 

where the upper sign pertains to the weakly compressible p < region and the lower 
one to the strongly compressible one p > |?, with f± giving two independent solutions 
of the eigenfunction equation (M — iuj)f = expressed by the Bessel functions: 



f±(p) = p v ^(^v^ p v (2 - 27) 

and with 

W = p a (h(p)d p f T (p) - f T (p)d P U(p)) (2.28) 

standing for their p-independent Wronskian. The eigenfunction / + ( /_) satisfies the sin- 
gular Dirichlet (Neumann) condition at the origin imposed by the limit when the trajectory 
noise is turned off for weak (strong) compressibility, see [19]. It follows that 

(e-*l {t<oo} ) = -d n[pi) (M - iu)-\ Po , Pl ) = (2.29) 

For uj — > the eigenfunction /_ reduces to a constant whereas /+ becomes proportional 
to p 1_a so that 

f 1 for p<4, 

^ = {(*)'" for P >i- ( " 0) 

Hence the exit time is almost surely finite in the weakly compressible regime whereas it 
is infinite with positive probability that depends only on ^ in the strongly compress- 
ible regime where the process p(t) is absorbed at the origin with the complementary 
probability. Such absorption corresponds to the coalescence of pairs of trajectories. We 
conclude that the Lagrangian flow is locally separating for p < d/^ 2 and locally trapping 
for p > dji 2 . 

For the conditional characteristic function, one obtains 





K 2 ^\ 


/fP^j 


J T b\ 




/fPi~ e ) 



= [yj — r-V—v- (2-31) 

The moments of the exit times may be derived from this expression by expanding the right 
hand side in powers of uj. In particular, one obtains for the conditional average of the exit 
time the result: 
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which reproduces in the £ — > 2 limit the n = 1 version of Eq. (2.24). The decay oc 
g-cCv^H) Q f the absolute value of the right hand side of Eq. (2.31) at large positive or 
negative uo guarantees that the exit time PDF Q(t; p , p) is smooth. Since it is zero 
for negative t, it must vanish with all derivatives at t — 0. More exactly, the decay 

~ e _6l v / H of the characteristic function (2.31) along the positive imaginary axis of uj, with 

h = 2(2-()- 1 (D' 1 )- 1 / 2 (pf^ )/2 - p { ^ 0/2 ), signals the behavior ~ e"« of Q(t;p ,pi) for 
small t. The analyticity properties of the right hand side of (2.31) imply the exponential 
decay ~ e~ b2t of Q(t; po, pi) for large t, with the rate b<i = i(2 — £) 2 D[pi~ 2 x 2 , where 
xq is the (real) zero of J^b(z) closest to the origin. In this respect, the exit time PDF 
Q(t;po,pi) behaves similarly for the weak and for the strong compressibility, the main 
difference between the two cases consisting in the missing probability in the latter case. 

For pi < po, the statistics of the time of exit through pi is related to the resolvent 
kernel of the generator Mp on the interval [pi,oo). For < £ < 2, the latter is given 
by a formula like (2.26) but with the overall minus sign, the cases po < p and po > p 
interchanged, and the functions f T , f± replaced by the Hankel functions 

/«(,) = P^^(^yi^) (2-33) 

for % — 1, 2, respectively. The square root in the argument of the Hankel functions should 
be taken with the positive imaginary part so that it is the eigenfunction which has a 
stretched exponential decay for large p. For the characteristic function of the exit time, 
Eq. (2.19) gives: 



c 



1 { ;<oo}> = TTrrf^. (2.34) 



, / (1) (Po) 
/ (1) (Pi) 

When uj — > 0, the eigenfunction becomes proportional to p 1_a if a > 1, i.e. if 
6<0 or p < 4=2 _)_ I anc i to a constant if a < 1, i.e. if b > or p > ^ + \. It follows 
that 

(W)-(ter 1 f "<g + f. (2.35) 
1 1 for p > ^ + I . 

For p < ^^ + | (i-e. for a > 1), the process p(t) remains forever in the interval (pi, 00) 

/ \ a— 1 , 9 . 

with probability 1 — (^-J > whereas for p > + ^ it exits through pi almost 
surely corresponding, respectively, to a locally transient and a locally recurrent Lagrangian 
flow. 

The conditional characteristic function is given by the expression 



<e") = r V , ' - (2-36) 
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Again, its absolute value decays as e for large |u;| implying that Q(t; po, pi) is 

smooth and vanishes with all derivatives at the origin. More exactly, the decay ~ e - b >-\^\ 
of (2.36) along the positive imaginary axis, where b\ is as for p 1 > p but with p and 

- b l 

pi interchanged, implies again the behavior ~ e ** of the exit time PDF Q(t; Po, Pi) for 
small t. Since H^\z) is a combination of z ±b with coefficients that are entire functions 
of z 2 (for non- integer b), it follows that ( e luJt ) c has the n th derivative over u at the origin 
if (and only if) n < \b\. That implies that for p\ < p the exit time PDF Q(t; po, Pi) has 
a power decay for large t, unlike for p± > p where it decayed exponentially: the large 
deviations of t are even more non-Gaussian in this case. 

For both pi > po and p\ < po, the Kraichnan model exit time PDF has for < £ < 2 
the scaling property: 

pQ{pt;p^p ,p^ Pl ) = Q{t;p , Pl ). (2.37) 

It follows that the moments of the exit time ( t n ) are proportional to p^ 2 ^ if ^ is kept 
constant. The same scaling implies the behavior (2.5) that characterizes what we have 
called the Richardson flows. This behavior may be also directly inferred from Eqs. (2.29) 
and (2.34). 
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FIG. 1. Phase diagram of the Lagrangian flow for the three-dimensional Kraichnan model. 

One may summarize the properties of the Lagrangian flow in the Kraichnan model in 
the phase diagram, drawn in Fig. 1 for three dimensions, with five phases that we list with 
their characteristics (assuming for £ > 2 finite L): 
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I. deterministic, Lyapunov, locally separating, locally transient, 
for £ > 2 and p < d/A ; 

II. deterministic, Lyapunov, locally trapping, locally recurrent, 
for £ > 2 and p > d/A ; 

III. stochastic, Richardson, locally separating, locally transient, 
for < £ < 2 and p < ^ + \ ; 

IV. stochastic, Richardson, locally separating, locally recurrent, 
for < i < 2 and ^ + \ < p < d/£ 2 ; 

V. deterministic, Richardson, locally trapping, locally recurrent, 
for < £ < 2 and d/£ 2 < p . 

These phases were essentially enumerated in [19] (with little stress put on the difference 
between phase III and IV). The characterization described above is closely related to 
the description of the phase diagram in [24], see also [30]. The notable difference is 
that, in order to characterize the phases in the non-Lipschitz case, refs. [24,30] used the 
dichotomic behaviors of the time of exit through pi in two limits: when p 1 — > with 
Po kept constant and when p — > with p 1 kept constant. Such behaviors enter the 
standard classification [15,6] of the one-dimensional diffusion p(t) on the half-line [0, oo[ 
with p = being, for £ > 2, a natural boundary and, for £ < 2, an entrance boundary in 
the weakly compressible phase III, a regular boundary in the intermediate phase IV and 
an exit boundary in the strongly compressible phase V. The use in the present paper of 
the small p behavior of the exit time PDF at fixed ^ in order to characterize the phases 
was motivated by the fact that such behaviors were both more amenable to analytic 
arguments in the presence of temporal correlations of velocities and more accessible to 
numerical simulations. 

As noticed in [30], see also [25], the solution for the PDF V(p ,p;t) corresponding 
to the singular Dirichlet boundary condition for M and coalescent trajectories, which 
pertains only to phase V if the flow is defined by adding and removing small noise, sets 
in already in region IV if we add no noise but first smoothen out the velocity fields at 
short distances and subsequently remove the smoothing. Physically, the first procedure 
corresponds to the vanishing Prandtl and the second one to the infinite Prandtl numbers. 
For 4^ + | < p < and well tuned Prandtl numbers, one may also obtain intermediate 
solutions that correspond to a "sticky" behavior of fluid particles [20]. The different 
limiting procedures give then rise to different boundary conditions that may be imposed 
on the generator M of Eq. (2.7) in the situation when p = is a regular boundary 
[15,6,31]. 
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3. GAUSSIAN VELOCITY ENSEMBLES WITH TEMPORAL CORRELATIONS 



The temporal decorrelation of the Kraichnan velocities is a simplifying feature that is 
quite unphysical since realistic turbulent velocities are correlated at different times. In the 
present paper we attempt to study the effect of temporal correlation of velocities on the 
behavior of the dispersion of a pair of particles in simplest ensembles of velocities with 
such correlations built in. More specifically, we shall consider the Gaussian ensembles of 
rf-dimensional velocities with mean zero and covariance 

(v%r) v'{f,r')) = D 2 J P»(k,p) *p. (3.1) 

There are three parameters in (3.1) not related to the choice of units: the spatial Holder 
exponent a, that we shall restrict to the interval (0, 1), the temporal exponent f3 taken 
positive, and the compressibility degree p G [0, 1] . Besides, there are three dimensionful 
parameters: D 2 of dimension length 2 ^ 1 " 01 ^ /time 2 , D s of dimension length 213 /time, and 
the integral length scale L. Similarly as in the Kraichnan ensemble, L may be taken 
to infinity for correlation functions of differences of velocities Av(t,p) whose statistics 
becomes scale invariant in this limit. The correlation time r c (p) of the velocity differences 
in ensembles given by Eq. (3.1) is equal to D^p 213 whereas the variance ((Av(t, p)) 2 ) = 
S(p) 2 behaves as D 2 p 2a . 

We shall be looking at the statistics of the 2-particle separation either using equation 
(1.2) with the ensemble (3.1) governing the Eulerian velocities, or using Eq. (1.4) with 
the ensemble (3.1) describing the quasi-Lagrangian velocities. It should be stressed that 
the two choices lead to two different models of Lagrangian flow. They exhibit different 
behaviors even for incompressible velocities where the equal-time velocity statistics of the 
Eulerian and quasi-Lagrangian velocities coincide. This should be contrasted with the 
situation in the Kraichnan model where the Eulerian and the quasi-Lagrangian velocities 
had the same all-time statistics so that it did not matter which one was modeled with 
the Gaussian ensemble (2.1). That the situation is different in the presence of temporal 
correlations is due to the manner in which sweeping by large scale eddies is taken into 
account. The 2-particle separation p involves only differences Av qL (t,p) of the quasi- 
Lagrangian velocities, see (1.4). The statistics of such differences has a regular L — > 
oo limit if we use the Gaussian ensemble (3.1) for the quasi-Lagrangian velocities. In 
this case the sweeping by scale L eddies does not effect the pair dispersion. On the 
other hand, p cannot be expressed in terms of the Eulerian velocity differences only, 
due to the dependence on the reference trajectory R(t), see (1.2). As a result, if we 
substitute the ensemble (3.1) for the Eulerian velocities, the dispersion statistics is still 
effected by the scale L sweeping and behaves in a singular way in the limit L — > oo. 
This singularity modifies also the short time behavior of the pair dispersion at fixed L 
in certain regimes, as we shall discuss below. The use of the synthetic ensemble (3.1) to 
describe turbulent velocities is in any case an approximation. It seems to render better the 
Lagrangian features of real turbulence if used to model the quasi-Lagrangian velocities. In 
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this case the large scale sweeping influences only the single particle statistics, but not the 
pair dispersion. We shall limit ourselves to this situation throughout most of this paper, 
dropping the subscript "gL" on the velocities. The exception is Sect. 7 where we discuss 
what happens when the Gaussian ensemble (3.1) is used to model the Eulerian velocities. 



FIG. 2. Phase diagram of the three different regimes of Lagrangian flow in the time-correlated 
velocities discussed here. The exponent a is the spatial Holder exponent and the exponent (3 
controls the behavior of correlation time vs scale. 

The first idea about the Lagrangian flow in the ensembles (3.1) of quasi-Lagrangian 
velocities may be gained by comparing the correlation time r c (p) to the eddy turnover 
time T e (p) = D 2 1 ^ 2 p 1 ^ a oc p/£(p). On the line a + 2/3 = 1 which, in particular, contains 
the Kolmogorov point a — (3 — | (see Fig. 2), both times have the same scale dependence. 
For a + 2(3 < 1 (domain C on Fig. 2), r c becomes much shorter than r e at large scales 
and much longer at small ones. We shall see that in this domain of parameters the pair 
dispersion should be described at long scales (for L = oo) by the Kraichnan model with 
rapidly decorrelating velocities and at short scales by a model with velocities independent 
of time (frozen). For a + 2(3 > 1 (domains A and B on Fig. 2), the relation between 
the correlation times is reversed and we could expect that the frozen model controls the 
large scale dispersion and the decorrelated one governs the short scale behavior. A similar 
picture underlined the phase diagram in a simple family of scale-invariant shear flows [26]. 
We shall study the asymptotics of the Lagrangian flow using scale transformations. Such 
transformations induce a flow in the plane of dimensional parameters of the model whose 
asymptotics is controlled by fixed points, as in the field theory renormalization group 
[1]. The perturbative renormalization group has been previously used in [2] to analyze 
the related scalar advection problem in the family (3.1) of velocity ensembles around 
{a, 13) = (—1,1) and there is some overlap of those results with our conclusions. The 
much heavier analysis of [2] was concentrated, however, on the aspects of advection related 








1 a 
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to finer details of the Lagrangian flow, see also [8]. Our point is that the analysis of the 
Lagrangian dispersion may be performed, at least to at certain extent, in a straightforward 
and nonperturbative way. 



The Kraichnan and the frozen model may be viewed as special limiting cases of the 
Gaussian velocity ensembles (3.1). The first one, with £ = 2 (a + (3), is obtained when 
D 2 , D 3 — > oo with ^ = D 1 kept constant as a consequence of the convergence 

D ^ k W e -\t\D 3 kf _^ 26{t y (32) 

L D3-+00 

At L = 00, existence of the limit for the correlation functions of Av requires that £ < 2, 
i.e. that a + ft < 1. Note that the convergence is fast at large wave number k, i.e. at 
small distances, and for long times, but it becomes slow for short times and, if L = 00, 
at large distances. The frozen model is obtained by taking D 3 — > with D 2 = const. In 
this case, 



-D 3 \t\kf 

but the convergence becomes slow at large k, i.e. at small distances, and for long times 



D 3 ^0 V ' 



Using the evolution equation (1.4) for the trajectory separation vector, we obtain for 
the mean rate of growth of the square of pair dispersion: 

j t ( P\t) ) = 2 j ( Av(t, p(t)) ■ Av(s, p(s)) ) ds . (3.4) 


Let us start by a naive mean-field-type approximate evaluation of the right hand side 
in the limit when p ~~ *■ 0- Such an evaluation should render correctly the behavior of 
{p 2 {t) ) in the stochastic regime. It is obtained by rewriting Eq. (3.4) as 

j t ( P 2 ) = 2T ((Avnt,p)), (3.5) 



where 



j(Av(t,p(t))-Av(s,p(s)))ds 

T = °- , , (3.6) 

7 (A^(t, P (t)) v 



has, if smaller than t, an interpretation of the correlation time of the Lagrangian velocity 
difference Av(t, p(t)). We may try to close Eq. (3.5) by assuming that T depends on 
the mean separation distance (p 2 ) 1 ^ 2 the same way as the correlations time r c (p) on p 
if T c (p) is smaller than t or as t otherwise: 

T w mm{0({ P Y), t} . (3.7) 
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Different domains in the space of parameters correspond to different choices of the minimal 
value (3.7) for T. As for the other term on the right hand side of (3.5), we shall again 
ignore the velocity dependence of p(t) putting 



(Avf) » 0((pY) (3.8) 

Using the above approximations for long times and L = oo, one obtains from (3.5): 

(p 2 {t)) « 0(t^) fora + 2/3>l (domains A and B) , (3.9) 

(p 2 (t)) « O(f^=0) for« + 2/3<l (domain C). (3.10) 

In the same way, we may estimate the short-time behavior in the limit po obtaining 

(p 2 (t)) « 0(*t===3) for a + /3<l<a + 2/3 (domain B) , (3.11) 

(p 2 (t)) « (9(tT=s) for a + 2/3<l (domain C) . (3.12) 

Note that, in agreement with (3.7), it is the smaller exponent that is chosen for large times 
and the bigger one for short times, a manifestation of a tendency of close trajectories to 
stay close. The region a + (3 > 1 (domain A in Fig. 2) which escapes the short-time 
estimates has been rigorously analyzed with the Gaussian ensemble (3.1) used to model 
the Eulerian velocities [14] and was conjectured to correspond to deterministic trajectories. 
We expect that also in the quasi-Lagrangian model the pair dispersion will concentrate in 
domain A at p = when p — > 0. This is consistent with the divergence of the predicted 
power in the short-time Richardson law (3.12) when a + (3 approaches 1 from below 
(i.e. from the domain B in Fig. 2). A similar divergence occurs in the weakly compressible 
Kraichnan model when we approach the Lipschitz regime £ > 2 from the non-Lipschitz 
one £ < 2. It signals there the passage from the power law to the exponential separation 
(2.10) of trajectories. 



4. SCALING ARGUMENTS 

The main aim of this note is to substantiate further the above conclusions based on the 
naive estimates (3.7) and (3.8). We shall also acquire an insight into the behavior of 
general moments of the pair dispersion and of the exit time and into the extent of the 
different Lagrangian flow regimes in the ensembles given by Eq. (3.1). In the study of the 
long- and short-time behavior of trajectories, it is convenient to consider their rescaled 
versions R^it) = p~ a R(pt) for appropriately chosen a. Since 

riff 

= v»(t,R»), (4.1) 

for the rescaled velocity v^(t,r) = p 1 ~ a v(pt ) pfr), the path R^(t) is a Lagrangian tra- 
jectory for Vp. There are special cases when the rescaled velocity differences have in the 
L — > oo limit the same distribution as the original ones for an appropriate (and unique) 
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choice of a. This happens for a = j^- both on the line a + 2(3 = 1 and in the frozen 



model and for a = 7^ in the Kraichnan ensemble. We infer that in those cases the pair 
dispersion PDF V and the exit time one Q are scale-invariant: 

/i CT P(/i>,//V;/it) = V( Po ,p;t), /1 /i>,/i>) = Q(t;p ,pi). (4.2) 

In particular, in the stochastic regime, the pair dispersion moments (p(t) n ), if finite in 
the po limit, behave as 0(t ncr ) for long times and become proportional to t ncr for all 
times when po —> 0. These conclusions fail in the deterministic regime, as we have already 
noticed in the Kraichnan model, see Eq. (2.15). In all regimes, the exit time moments 
(t n 1{<<oo} ) (if finite) are proportional to p^ CT if ^ = 7 is kept constant. 

Out of the line a + 2(3 — 1, the scale invariance of the Lagrangian dispersion is broken 
but in a predictable way, as we shall see. The crucial observation is that the rescaled 
velocities v^{t,r) are distributed with the 2-point function (3.1) with Di replaced by 
D i fl and L by L M , where 

D 2 ,» = p 2 ^ {1 - a ^D 2 D 3t , = p 1 - 2 ^ D 3 , L, = p~°L . (4.3) 

For p tending to infinity or to zero (i.e. when exploring the long-time or the short-time 
behavior of the particle separation), we may choose o so that the distribution of the 
rescaled velocity differences at L = 00 tends to the Kraichnan or to the frozen model ones 
(with the notable exception of the p — * limit in domain A). 

Consider first a — > 00 at L = 00. Taking a = 7m — ^ — ^ fixes the ratio 2 ,P 2 ' M with 
D i41 ^00 if a + (3 > 1 (domain A in Fig. 2) or if a + 2/3 < 1 (domain C in Fig. 2). 
The latter case leads to a non-singular Kraichnan ensemble of velocity differences with 
£ = 2(a + (3) whereas the former one does not (it would correspond to £ > 2, L — 00). 
We may then expect that 



Jim p CT P(p CT p ,p CT p;pt) = V Kr (p ,p;t) 
lim p Q(pt;p a p ,p' 7 p 1 ) = Q Kr (t; p , pi) 



} for (J =W^F) and « + 2/ 3<1 ( 4 - 4 ) 



where V Kr and Q Kr pertain to the Kraichnan model with £ = 2(a + (3) . This is indeed 
consistent with the scaling properties of the Kraichnan model dispersion. 

Taking a = fixes D 2 ^ with D 3tfl — > if a + 2(3 > 1 (domains A and B in Fig. 2). 
We then expect that 



lim p°V(p°p ,p a p;pt) = V fr (po,p;t) 



[1— >oo 



lim p Q(pt;p a p ,p a pi) = Q h {t;p ,pi) 



J for a = and a + 2(3 > 1 , (4.5) 



where V fr and Q fr stand for the PDF's of the frozen velocity model with Holder exponent 
a. Note again the consistency with the scaling properties of the Lagrangian dispersion in 
the frozen model. 
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Inquiring about the short-time asymptotics of the trajectory dispersion reverses the 
asymptotics. We should then have 



^° \ for a = j-^ and a + 2(3 < 1 (4.6) 

lim p Q(pt; p° p , p° p x ) = Q /r (t;p ,Pi) J 

(i.e. in domain C in Fig. 2) with the same value of the Holder exponent a, and 

lim ^V.,/^) =V^,p;t) , fma= md (a + P<l . 

lim i/, Q(pt;v"m,»°Pi) = Q ,l> (*;Po,Pi) J 2<1 "°"" ) l« + 2/3>l 

(i.e. in domain B in Fig. 2) with £ = 2(a + /3) for the Kraichnan model. Again, this is 
consistent with the scaling of the limiting PDF's. In summary, the scale invariance of 
the statistics of the pair dispersion and of the exit time, although broken away from the 
a + 2/? = 1 line, should be restored at long and short times. 

It has to be stressed that the relations (4.4) to (4.7) are conjectural. The PDF's V and 
Q are complicated nonlinear functionals of the quasi-Lagrangian velocity statistics and 
the conjectured relations assume their continuity in an appropriate topology, which is not 
obvious. In particular, since the convergence of the rescaled velocity covariances to the one 
of the Kraichnan model is very slow at long distances, there is a potential threat for the 
corresponding convergence of the rescaled PDF's V(po,p;t) and Q(t; p , pi) with pi < p 
coming from the contribution of trajectories that venture far apart, if such contributions 
are important. Similarly, the slow convergence to the frozen model at short distances 
could create problems for the corresponding convergence of the rescaled pair dispersion 
and exit time PDF's, for the latter if p\ > p . Whether such effects invalidate some of the 
conclusions (4.4) to (4.7) could be, in principle, studied in perturbation theory around the 
Kraichnan or frozen model. The question of the convergence of rescaled dispersion PDF 
to the Kraichan model has been recently rigorously studied in [12,13]. We shall further 
discuss the non-uniformity of the convergence of the exit time PDF to that of the frozen 
model in Sect. 6. 



An important question concerns the phase diagram of the Lagrangian flow for the 
model (3.1) of quasi-Lagrangian velocities. As mentioned above, domain A is expected to 
correspond to the deterministic Lagrangian flow. The rate of separation of close trajecto- 
ries in this domain (Lyapunov or Richardson flow? locally separating or trapping? locally 
recurrent or transient?) is also an open problem. Inside domains B and C, we may try 
to use the conjectured convergence (4.4) to (4.7) of the rescaled pair dispersion and time 
exit PDF's to characterize the nature of the Lagrangian flow. This will require even more 
care since some uniformity of the limits will be needed. 

First, we may argue that, inside domains B and C, weak compressibility p < 4(a ^ 2 
implies that the Lagrangian flow is stochastic. The argument assumes that the dichotomy 
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"deterministic versus stochastic" may be still characterized as in Sect. 2, i.e. by the behav- 
ior of V(po,p;t) in the limit p — > 0. It goes as follows. Suppose that the relation (2.2) 
holds at some point inside B or C (for all t). Then, obviously, also 

lim p°V(p a p ,p a p;pt) = 5{p) (4.8) 

po->0 

for all p. We expect the convergence (4.4) in domain C and (4.7) in domain B, both 
resulting in the Kraichnan model PDF's, to be uniform at short distances and hence to 
commute with the po limit. We may then infer from Eq. (4.8) that (2.2) holds also 
for the limiting Kraichnan PDF so that p > d/C, 2 where £ = 2 (a + /3), which implies the 
assertion. 

The analogous argument may be applied when 

lim / Q(t;p ,7p )dt = lim / p Q(pt; p u p , p a 'yp ) dt = 1 (4.9) 

leading to the predictions that inside the domains B and C, the Lagrangian flow is locally 
trapping if p > 4 ( a +p\2 and locally transient if p < 4 ^~^ + \- The first claim seems 
somewhat more trustable since it does not involve large separations where the convergence 
to the Kraichnan ensemble is slowed down. Both require additionally that no mass in the 
exit time PDF escapes to infinity during the p-limits (4.4) and (4.7). This should not pose 
a problem since convergence to the Kraichnan model becomes very fast at long times. 

In the stochastic regime, the convergence (4.4) and (4.7) should imply for the (positive) 
pair dispersion moments the behavior 

n 

(p n (t)) = 0(tw-°-n) (4.10) 

for long times in domain C and for p — > and short times in domain B, in agreement 
with the naive mean-field results (3.10) and (3.11). As for the behavior of the moments 
of the exit time through 7p , the relations (4.7) and (4.4) should imply that 

(i n W) = O(p 2 n(1 - a -V) (4.11) 

in domain B for small po and in domain C for large ones, irrespectively of the character 
of the Lagrangian flow. Similarly, we should obtain in domain B the convergence 

lim J f(t) Q(t;p°p ,p° 1Po )dt = /(0) J Q Kr (t;p , Pl )dt (4.12) 
characteristic of the Richardson flow in our terminology. 

Similar use of the conjectured convergences (4.5) and (4.6) to the frozen model PDF's 
in order to argue about the Lagrangian flow dichotomies and the evolution of dispersion 
and exit times moments poses two problems. First is the poor knowledge of the flow 
behavior in the frozen model, see, however, Sect. 5. Second, even more serious one, is 
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the non-uniform nature of the convergence that becomes slow at small distances and long 
times. It is still plausible, however, that the convergence (4.5) is uniform enough as to 
imply that 

(p n (t)) = 0(tT&t) (4.13) 

for long times and for sufficiently high n in the stochastic regime in domains A and B. 
Similarly, relations (4.5) and (4.6) may still imply that for sufficiently negative moments, 

(t n ) = O{p n {l - a) ) (4.14) 

in domains A and B for large po and in domain C for small ones, see Sect. 6. The same 
way, the convergence (4.12) should still hold in domain C for test functions / decaying 
fast at infinity and Q Kr replaced by Q fr . 



Summarizing, we predict, with various level of confidence, that the Lagrangian flow is 



deterministic 
in domain A 



stochastic, Richardson, locally transient 
in domains B and C for p < 4 f~+ 0) + \ , 

stochastic, Richardson 
in domains B and C for J~+ 3) + ^ < p < 



d 



4(a+/3) 2 4(Q+/3) 2 ' 



Richardson, locally trapping 
in domains B and C for p > 



4(a+/3) 2 ' 



The degree of confidence of the predictions depends on which of relations (4.4) to (4.7) 
was used in the argument and with what uniformity assumptions. The predictions are 
consistent with the intuition that increase of the compressibility degree p enhances the 
trapping of the fluid particles. It will be interesting to confirm (or infirm) them by further 
analytic arguments and by numerical simulations. Note that, in particular, we expect 
that in domains B and C the incompressible Lagrangian flow is stochastic, Richardson and 
locally transient, and that, if the dimension d > 4, it stays such, whatever compressibility. 



5. ONE-DIMENSIONAL FROZEN ENSEMBLE 



The frozen model was left out from the discussion of the phase diagram in the last section. 
Our arguments were based mainly on the convergence of the rescaled velocity ensemble 
(3.1) to the Kraichnan model and such convergence is, of course, absent for the frozen 
model. One may expect appearance of discontinuity in the character of the Lagrangian 
flow in the frozen model limit D 3 — > which is is very non-uniform at short distances and 
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long times, leading to a strong enhancement of trapping. This effect will be analyzed in 
next Section. Here we shall try to find out what happens in the frozen velocities which, in 
general, are the hard to analyze. One case where some analytic results may be obtained is 
the one-dimensional Gaussian model with the Holder exponent a — \ and 2-point function 

/, p ik(x—y) Jl, 
= i^-,*-,,/^ (51) 

where, for simplicity, we have set D 2 = 1 (what may be always achieved by rescaling 
v I— > y/D^v). A simplifying feature of this case, studied already in [29], is that v(x) forms 
a stationary Markov (Ornstein-Uhlenbeck) process with the generator 

£=~ 1 2& + ii (5-2) 
corresponding to the Focker-Planck harmonic oscillator Hamiltonian 

H = e -£e- = --— + _--. (5.3) 



The velocity v(x) with fixed x is distributed according to the invariant measure of the 
process 



,,2 



du(v) = -= e l dv . (5.4) 

\/ttL 



The transition probabilities of the process are 



p(v ,v;t)dv = e- tc (v ,v)dv = 1 e Hi-^) dv . (5.5) 

y iTL(l-e^ 2t ' L ) 

In the limit L — > 00, the velocity difference Av(x) = v(x) — v(0) becomes the one- 
dimensional two-sided Brownian motion w(x). The quasi-Lagrangian Eq. (1.4) for the 
trajectory separation takes then the form of the steepest descent equation 

I = «(P) = (5.6) 

x 

in the potential (f>(x) = — J w(y)dy. The solution p(t) slides to the bottom of the 



potential well in which the initial point p = p(0) is situated, i.e. to the closest zero p + 
of w(p) to the right of p if w(p ) > or p_ to the left of p is w(po) < 0. The only 
difference with the case of smooth potential with wells approximately quadratic around 
typical minima, is that, as we show below, the solution will arrive to the bottom of the 
well in a finite rather than infinite time. This is due to the roughness of the Brownian 
motion. After arriving at the bottom, the solution will stay locked there in subsequent 
times. 
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We may restrict ourselves to the case p > 0. Suppose also that w(p ) > 0. The first 
value p + to the right of p such that w(p+) = is finite with probability one. The time 
t + that the solution p(t) of Eq. (5.6) starting at p at time zero takes to reach p + is 

U = (5.7) 

Let us compute the expectation (t + l{iu(p„)>o,p + <p 2 }) °f times t + over the Brownian paths 
w(p) such that w(p ) > and p + < p 2 for certain p 2 > po- It is equal to 

P2 00 2 °° 9 OO „ 

dp e 2p o ; / e 2 (p-po) — e 2 (p-po) — / e 2 (p 2 -p) — . (5.8) 

7 y V^,/ V / y /2n(p-p )w 7 \f2-wip2-p) 

Po 

The origin of this formula is straightforward. The probability that w(po) belongs to 

w 

[w ,w + dwo] is e The one that w(p) belongs to [w,w + dw] without passing 

through zero between p and p is e 2 (p-po) — e 2 (p-po) ) , , given that u>(p ) — 

V / v /27r(p-po) 

w > (it is expressed by the heat kernel with the Dirichlet condition at w — 0). Finally the 
last integral on the right hand side of (5.8) is equal to the probability that the Brownian 
trajectory crosses zero between p and p 2 , given that w(p) = w. In Appendix 1 we show 
that 

(t+Mw(po)>o,p + < P 2}) < C+ pl /2 In , (5.9) 

where C + is a dimensionless constant. This proves that the time t + is almost surely finite 
although the unrestricted mean ( t + l{«,(p )>o} ), given by the p 2 ^ oo limit of (5.8) under 
which the last integral on the right hand side tends to one, diverges. The divergence is 
due to the contribution of the Brownian paths that travel far before falling back to zero. 
In the Ornstein-Uhlenbeck process with L < oo, the weight of such paths is suppressed 
and it is not difficult to show that (t+ l{ w {p )>o} ) is finite then. 

Similarly, let w(po) < and < p_ < po be the first value to the left of po such 
that w(p-) = 0. The time i_ that the solution p(t) of Eq. (5.6) starting at po at time 
zero takes to reach p_ is given by Eq. (5.7) with t+ replaced by t- and p + by p_. The 
expectation value {t- l{ w (p )<o} ) is given by 

po o o „ „ 

— / dp I e 2 p 7 = / e 2 (po-p) - e 2 (po-p) (5.10) 

J J w J \ J x / 2 7r(po-p) 



-oo 



which is easily seen to be finite, e.g. by bounding the last integral by ^ n (po-p) \ w \- ^ e 
infer that 

( t- l{ w (p )<o} ) = C-Po /2 , (5.11) 
where C_ is another dimensionless constant. 
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A. Exit time 



As for the exit time that the process p(t) takes to grow from p to p 1 > p , it is finite if 
and only if p(t) is not stuck at a zero p± < p x of w(p), i.e. if w is positive on the interval 
[po, Pi)- It is then given by the formula 



Po 

The probability that t < oo is given by 



w{p) 



OO 2 00 o 2 

(l {t<oo} ) = Je «»-p=J (e ^> -e ^-,o»)^_. (5.13) 
00 v 

Note that the last expression depends only on ^, is smaller than 1/2 and tends to zero 
when po with px fixed. With the complementary probability the solution p(t) 
starting from p > will never reach p 1 , i.e. the exit time t is infinite. The averages of 
positive powers of the exit time constraint to be finite are expressed by the relation 



(t n l {t <oo}> = n\ J dp'---dpW Je~% 



2 

dwo 



Po<p'<---<p (n) <pi 

7J / ( e UpW-pV-V) - e 2(pW-p(*-D) V 



i=i 



du>(0 



/ 



)2 („W +roi ) 2 

f e 2(pi-p (n) ) - e -,(")) ) — dwi (5.14) 

V / v /27r(pi-p(")) 

with p(°) = po an d w< ' ' ) = Wo- It is easy to show t hat the ex pression on the right hand 
side is finite. Indeed, bounding the last integral by \f~^^^) w<yn) as i n estimating (5.10) 
and proceeding further the same way, we obtain the inequality 

nil r , % 



V(p'-Po)-(Pi-P (n) ) 
Po<p'<--.<p (n) <pi 



n-1 



where the last line results from the inductive calculation of the p^ integrals. As a conse- 
quence of the scaling properties of the right hand side of Eq. (5.14), 

(t n l {t< oo } ) = C n pT = C n r n e (5.16) 

where C n are dimensionless constants depending on the ratio 7 = ^ > 1 and r e = 
(p /D 2 ) 1 ^ 2 is the eddy turnover time at scale p (recalled that we have set D 2 — 1 
above). 
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The bounds (5.15), which respect the above scaling, imply that the characteristic 
function (e lujt l{t <0 o} ) is entire in u and that the large t decay of the PDF Q(t; p , pi) 
is at least oc e~ olyt \ Due to Eq. (5.12), the characteristic function may be expressed by 
the path integral 

(^W) = bfe-rtj^fero V w (5.17) 

o 

over the paths [po, Pi] 3 p i— > w(p) G (0, oo) such that w(p ) = w , with Af being 
an appropriate normalization factor. The expression permits to evaluate the large 
behavior along the positive imaginary axis of uj by the semi-classical calculation. The 
extremal trajectory w(p) describes a motion of a unit mass particle climbing up in the 
potential — ^ until a total stop. It satisfies the equation 

sin i^o . -\ Pi 



(<p + sin ip cos (p) = (5.18) 



COS 3 po v ' ' ' ' PO PO 



for cos 2 y? = ^^y, cos 2 (p = ^j^y and (p between zero and (p < § • In particular, tp 
is determined by Eq. (5.18) with ip — ip and p = po- It depends only on the ratio ^. 
The action of the classical trajectory is S = \uj\ 2 / 3 pl^ 3 s with a dimensionless constant 
_ 3 *f?/ s . inl/3 y " growing with the ex . The decay ~ e -50 of the characteristic function 

u 2 1 /' 5 cos <po ° ° po 

4 3 . 2 

along the positive imaginary axis corresponds to the small t behavior ~ e ~27 poS ° = 
e~27 s o(W*) of the exit time PDF Q(t; p , pi) (up to powers of t). One could expect that 
for short times the trajectories move almost ballistically: 

pit) ~ p + ^(po)t (5.19) 

which would give t ~ , p ? and a short time tale ~ e 2p o of the exit time PDF. 
This reproduces well the power of t in the exponential but not the coefficient. The latter, 
divided by po, agrees only to the order (^ — l) 2 with the correct expression t^Sq. 

For uj on the imaginary axis, the path-integral on the right hand side of Eq. (5.17) may 
be re-expressed in the operator language via the Feynman-Kac formula. Using also the 
invariance of the Brownian motion under the scale transformations w(x) i— > |c<j| _1 w(|a;| 2 x), 
we obtain the identity 

OO 2 oo 

<e^l {t<oo} > = /e-^-=^-/e-H 2 ^-^ ± ( Wo ,^)^i, (5.20) 



\J2irpa M 




where the operator 
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on the interval [0, oo) is a 1-dimensional Schrodinger operator with Dirichlet boundary 
condition at zero. The signs pertain to the positive or negative imaginary cj-axis and result 
in a repulsive or an attractive potential, respectively. An explicit expression for the kernel 
of the exponential of /C_, see Appendix 2, shows the growth of the characteristic function 
for uj = —i\u)\ dominated by the lowest bound state contribution ~ e~'^' {pi-po)e 
E = -\. This implies the decay ~ e~T^) t2 = e -^r (^(po)) 2 of the pDF Q(t; Po , Pl ) 
for large t. 

The above analytic predictions may be used to validate the numerical method to be 
applied later for the cases a 7^ 1/2, where such rigorous results are not available. The 
velocity field is generated by using a straightforward Fourier method, i.e. generating 
the Fourier modes by a standard Gaussian random number generator and transforming 
back to real space by Fast Fourier Transform. The resolution which we used was 2 20 . 
The initial separation is set at 20, 000 and we measure the PDF of the time taken to 
reach 5 times the initial separation. The curve is shown in Fig. 3. The agreement with the 
superposed predictions indicates that the choice of the resolution and the initial separation 
are appropriate to avoid contamination by periodicity and/or discretization effects. 
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FIG. 3. The PDF of the exit times for the frozen flow with a = 1/2. The dotted curves are 
the analytical predictions at small and large times, respectively. 



The time in which p(t) decreases from po 1° a positive value p\ < po is finite and 

still determined by Eq. (5.12) if and only if w < on the interval (pi,p ]. It is easy to 
see that 

Q{t; Po , Pl ) = Q(t;pi,p ) (5.22) 
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since for pi < p the characteristic function {e tuJt l{t <0 o} ) is given by the expression (5.17) 
with p and p 1 interchanged. 

The scaling property (4.2) with o = 2 and the fact that the exit time is infinite with 
a positive probability depending on ^ imply that 

lim Q(t;p 0)7 p ) = c ( 7 )5(t) (5.23) 

po— >0 

for any positive 7 ^ 1 with < 0(7) < 1. It follows that the Lagrangian flow is 
Richardson, locally trapping and locally transient in the terminology of Sect. 2A. 



B. Pair dispersion 

The pair dispersion in the model is closely related to the exit time. The reason is that the 
solution p(t) of Eq. (5.6) never changes direction. As a result, the solution that starts at 
Po at time zero satisfies p{t) > p\ > p or p(t) < p\ < p if and only if it exits through 
pi in time shorter than t. It follows that 

00 t 

J V(p ,p;t)dp = J Q(s; p , p ± ) ds for p 1 > p , (5.24) 

pi 
pi t 

J V(p , p;t) dp = J Q(s; p , pi) ds for p 1 < p (5.25) 


or that 

t 

V(p ,p;t) = Td p J Q(s;p ,p)ds for p < p . (5.26) 


00 

The right hand side of Eq. (5.24) is bounded by / Q(s; po, pi) ds which tends to zero when 



Po — > for pi > fixed, see Eq. (5.13). It follows that 

hm V(p ,p;t) = 5(p), (5.27) 

in agreement with the fact that the trajectories are determined by the initial condition. 
Hence the Lagrangian flow is deterministic. Note the important difference with the phase 
V of the Kraichnan model, see Fig. 1. There the trapping effects resulted in the coales- 
cence of fluid particles signaled by the singular form (2.13) of the pair dispersion PDF. 

With the use of symmetry (5.22), the left hand side of Eq. (5.25) may be bounded by 

00 

/ Q(s] Pi,Po) ds and hence tends to zero when pi — > for fixed po > 0. As a result, the 


PDF P(po, p; t) cannot have a contribution proportional to <5(p) for p > 0. It is, instead, 
regular in p. The reason is that in the quasi-Lagrangian model of particle separation 
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(5.6), the solutions pit) are trapped with probability 1 at positive zeros of w (there is 
an infinity of such zeros) and never arrive at p = 0. This seems to be an artifact of the 
frozen one-dimensional model but it serves as a warning that the behavior of trajectories 
in the time-correlated velocities may be richer than what was observed for the Kraichnan 
model, with a possible occurrence of phase 

VI. deterministic, Richardson, locally trapping, locally transient 

characterized by a combination of properties that did not occur in the time-decorrelated 
model. 



C. Case with ol^\ 

Several of the features analytically established for the one-dimensional frozen model with 
the Holder exponent a — \ may be generalized to the case of general a in the inter- 
val (0,1), although the absence of the Markov property in the process v(x) makes the 
arguments more difficult, see [21]. In the limit L — > oo, the equation for the trajectory 
separation still has the form (5.6) with w(x), x > 0, being (upon a right choice of the 
normalization constant D 2 ) the two-sided fractional Brownian motion (fBm), i.e. the 
Gaussian process with mean zero and 2-point function 

(w(x)w(y)) = \[x 2a + y 2a - \x - y\ 2a ) = G(x,y) (5.28) 

for x, y > 0. The 2-point function (5.28) is a kernel of a positive operator on the half-axis 
[0, oo) that we shall denote by G. Note the scale invariance under w(x) h- > p~ a w(px) of 
the fBm. The basic result of [27] asserts that the probability that w(x) < w for w > 
and all x in the interval (0, po) behaves like O((pow 1,/a ) 1_a ) for large values of pow 1 ^ a . 
Similarly as for the Brownian motion, the fBm lives on continuous trajectories and has 
zeros in any interval (0, po) or (po,oo). One may again prove that, with probability 1, 
the solution p(t) of Eq. (5.6) starting at any p > arrives in finite time at the closest 
zero p± to the right or left of p , see [21]. 

As in the case a — \, the exit time t through p 1 > p is finite if and only if w > 
on the interval [p ,pi) and the probability ( l{t<oo} ) of such an event depends only on 
^ and tends to zero when po — > 0. It should be again possible to extract the behaviors 
of the exit time PDF <2(£;po,Pi) for large and small time by looking at the large |u;| 
behavior of 

pi 

<e ±M % <oo} > = (e "o l {w>0on[fKhPl)} ), (5.29) 
where the last expectation is with respect to the Gaussian measure of the fBm w. For the 
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negative sign, this expression should be still dominated for large |u;| by the semi-classical 
contribution ~ e~ So . The classical trajectory w(p) = \uj\ 1 / 3 p 1+2a ^ 3 (Gu)(-^), where u is 



a function that does not vanish only on the interval (1, ^) and such that 

u = (Gu)- 2 (5.30) 



there. Note that it follows that w(p) > for p > since G(x,y) > except for x, y — 0. 
The action of the classical trajectory is S = \uj\ 2 ^ 3 p^ 1 ' ^ 3 s for 



pi/po 

so = 



J [jG^o) + \u(x) (Gu)(x)] dx. (5.31) 



4 2(1 — a) 3 ,_2 

Such a semi-classical dominance implies again the small time tail ~ e~^ p ° s ° = 
e -o((Te/t) ) Q f ex j|- ^j me PDF Q(t; po, pi) (with the eddy turnover time r e = 

D 2 1/2 pl a ). 

With the use of the scale invariance of the fBm, one may also absorb the \u (-dependence 
of the characteristic function into the length of the p-interval in Eq. (5.29): 

pi 

(e ±|w|t l {t <oo } ) = (e 1 {W>0 o„ [P>[)} ) , (5.32) 

where p\ = \uj\ rr * pi. For the positive sign in Eq. (5.32) and a — |, we have used 
the Feynman-Kac formula in order to extract the extensive behavior of the right hand 
side for large |u;|. For other values of a E (0, 1) such a formula is not available but the 
expectation (5.32) may be viewed as the partition function of a one- dimensional continuous 
spin system with long-range 2-spin interactions decaying as distance~ 2( - 1+a ^ and with 
partially confining (as opposed to the case with negative sign) single-spin potential. It 

is plausible that the extensive behavior of the partition function ~ e~' w ' a (pi-pa)E a f or 
large \uo\ with free energy density Eq < persists for such systems. Such a behavior 
would result in the long-time tail ~ e y \(pi-po) E o\' ' = e~ U{{ /Te> > of the PDF 
Q(t; po, pi). The validity of this prediction is confirmed by Figs. 4 and 5 where we present 
the PDF's of the exit times for the two cases a = 0.4 and 0.75. The numerical simulations 
are realized by the same method previously validated in the case a = 1/2. 



The statistics of the time of exit through p\ < po is again obtained from that for p\ > 
Po by interchanging po and p\ and the pair dispersion PDF is still given by Eq. (5.26). 
The remarks about the phase VI type behavior of Lagrangian flow carry over from the 
case a — |. 
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FIG. 4. The PDF of the exit times for the Id frozen velocities with a = 0.4. The dotted curve 
is a fit of the form e - const - . 




0.2 0.4 0.6 0.8 
times 



FIG. 5. The PDF of the exit times for the Id frozen velocities with a = 3/4. The solid curve 
is a fit of the form e ~ const - . Since the curve gives the visual impression of an exponential 
decay, we also plot in the inset In | ln(Q(t; po, pi))| vs logt and compare the slopes 0.75 and 1 to 
validate the former. 



30 



6. EFFECT OF LONG TIME VELOCITY CORRELATIONS 



The presence of permanent trapping of trajectories in one-dimensional frozen ensemble, 
leading to events with infinite exit time through pi = 7p , should also occur in higher- 
dimensional frozen ensembles obtained by setting D% = in Eq. (3.1). It is plausible that, 
at least in the presence of compressibility, the sets of velocities trapping trajectories in a 
given region have positive probability. If, however, one reintroduces finite temporal corre- 
lations of velocities by taking positive D 3 in (3.1), the particles will eventually be released 
from traps after the time evolution changes the velocity field configuration, i.e. after time 
of order r c (p) = D^p 2 ^ for traps of size p. If for D% > the exit times are finite so that 

oo oo 

/ Q(t) dt — 1, then the missing probability 1 — / Q(t) dt > for D 3 = should be built 
o o 

from the long-time tails of Q at positive D 3 , as the latter is taken to zero. We may then 
expect that, loosely speaking, for fixed po, 7 and _D 2 , 

Q(t;p ,7Po) -T-iQ^/rcT) for t/r c > 1 . (6.1) 

with r c = r c (po) and the function Q°° having a (Poissonian) exponential tail. Let us 
concentrate on the one-dimensional situation where the frozen exit-time PDF is expected 
to have a stretched exponential tail ~ q- ^/^) 1 *") with the scale set by the eddy turnover 
time r e = D 2 1 ^pl~ a . Note that the latter attains the value of order r c _1 for t — r' with 
r'/r e = (9((ln(r c /r e )) a ) = 0(\ lnD 3 \ a ). The "minimal" scenario would be that the frozen 
PDF passes into the form of (6.1) around t = t'. More precisely, we may postulate the 
convergence 

Q(t;po,ipo) l{t<r'} — > Q fr {t;po,ipo) 

L>Z— »U 

(6.2) 

r c Q(sr c ; po, 7Po) 1{ s >t'/t c } — > Q°°{s;i) 

in a strong enough sense. The above relations imply that the large exit time behavior 
becomes self-similar for small D 3 , with the characteristic scale equal to r c , with no in- 
termediate regime between the frozen type behavior and the self-similar tail. It is also 
possible that a different intermediate regime sets in between times of order r' and r c , 
with t' depending differently on _D 3 . 

The minimal scenario would make explicit the large t non- uniformity of the conjectured 
convergence of the exit time PDF's, see (4.5) and (4.6). Recall that the conjecture was 
based on the scaling relation that may be rewritten as the identity 

P - a Q(p^;Po,7Po)L 3 = Q(t;l,7)L 3(P0) (6.3) 

for D 3 (p ) = pl~ a ~ 2l3 D 3 and D 2 unchanged, see (4.3). When p — > oo in domain A and 
B and po — > in domain C then Ds(pq) — > so that we fall into the situation considered 
in the scenario (6.2). The latter becomes then the assertion that 



31 



Pl a Q(pI% po, 7Po) i{t<r'(po)} — >, 2 /r (*; 1, 7) 

Po-> n 



(6.4) 



A$ V<T 2(^3 VIT; Po, 7Po) l{ s >D 3 (po)r'( P( ,)} — : Q°°(s; 7) 



po- 



for r'(p ) = O(|lnp | a ). 

One of the consequences of such a limiting behavior would be a bifractal scaling with 
Po — ► oo or po — > of the moments of the exit time. Indeed, 

oo t'(po) 

jt n Q(t;p , 1Po )dt = pt a)n J t n f%- tt QU- tt t;p ,'yp )dt 
o o 

oo 

+ D- n pf n J s n D~ 1 pf Q(-D 3 1 pfs; po, 7Po) rfs . (6.5) 

D[i(pa)r' (po) 

In the regime of extreme values of po, the first term on the right hand side behaves 
like p ( o >n ft n Qf r (t;l,7)dt if we assume (6.4) with a sufficiently strong convergence. 



o 

oo 



Similarly the second term would behave like -Dg n p / s n Q°°(s; 7) ds. If Q°°(s, 7) is 

D3(po)t'(po) 

integrable at zero, the first term dominates for negative n and the second one for positive 
n when it behaves as C(po )■ The "minimal" scenario (6.2) would then imply that 

'O(p n {1 - a) ) for n<0, 

(t n ) = (6.6) 
.^(Po^) for n>0 

for large p in domain A and B and for small p in domain C. Such a bifractal behavior 

is, of course, consistent with the earlier conjecture (4.14). The missing probability in the 

00 

frozen case would be given by / Q OD (s; / -f) ds. Even in the presence of an intermediate 



regime in the exit time PDF, the scaling (6.6) should set in for |n| 3> I. The convexity 
(concavity) of the large (small) p exponent as a function of n, together with its vanishing 
at n — 0, would then impose the behavior (6.6) for all n. 

We have tested the scenario (6.2) numerically. The one-dimensional velocity field 
with temporal correlations was obtained again using the Fourier method. Each Fourier 
mode Vk(t) was generated by integrating the corresponding Uhlenbeck-Ornstein differential 
equation: 

where r(k) scales as k~ 2/3 , the energy spectrum of the field is denoted by E{k) oc A; -1-2 " 
and dW(t) is a standard Brownian motion. The stochastic differential equations (6.7) 
were integrated by using a simple Euler scheme of order | [22]. 
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FIG. 6. The PDF of the exit times for the flow with a = 0.75 and (3 = 0.3 at six initial 
separations p = 400, 800, 1200, 1500, 1800, 2100. The resolution is 32768. 
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FIG. 8. The same PDF as in the previous figure but plotted now as f?§ Q, vs p^ 2 ^t to display 
the collapse at large exit times. 



The parameters of the flow where a = 0.75 and (3 = 0.3. The exit times for six different 
initial separations was measured. Their PDF's Q(t; p , 1.05p ) are shown in Fig. 6. Ac- 
cording to (6.4), by plotting p\~ a Q versus Po~ 1 t all the curves should collapse at small 
exit times, as confirmed in Fig. 7. Furthermore, the prediction (6.4) at large exit times 
is verified in Fig. 8 by collapsing the long-time parts of the curves by plotting p Q ver- 
sus Pq 2/3 1 Note that the curves for the two smallest values of po are not collapsing, in 
agreement with the previous arguments predicting that the asymptotic behavior sets in 
at large p ' s - The simulations are consistent with the absence of an intermediate regime 
but do not really allow to exclude such a possibility. 



7. EULERIAN GAUSSIAN VELOCITIES: SWEEPING EFFECTS 

Let us discuss how the Lagrangian flow changes if the Gaussian velocity ensemble 3.1 is 
used to model the Eulerian velocities rather than the quasi-Lagrangian ones. As already 
mentioned, the main difference is that, unlike for the Kraichnan model, the separation of 
two trajectories is not any more a function of velocity differences only and it is strongly 
influenced by large-scale eddies or the so called sweeping. This effect grows with growing 
integral scale L and we shall attempt to study its large L asymptotics. It seems to be 
stronger for small values of (3, i.e. for velocities that are almost frozen at small distances. 

The r.m.s. value of velocity in the ensemble (3.1) is proportional to L a , i.e. it becomes 
large for large L. On the other hand, the r.m.s. equal-time velocity differences on scales 
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much smaller than L are of the order distance . In particular, on the scales ~ L a they 
are of the order L a <C L a . Rewriting the trajectory equation (1.1) as 2 

t t 

R(t) = Jv(s,0)ds + J[v(s,R(s))-v(s,0)]ds, (7.1) 



we may expect that, for fixed t, the first integral is of the order L a and the second of the 

2 

order L a C i a . For bounded times, the first integral should then give the term of the 
order L a of the solution and the second one, with R(s) replaced by the approximation 
~ L a , the term of the order L a . More precisely, let us observe that the Gaussian process 
with the components L~ a v(t,0) and L~ a [v(t,L a r) — v(t,0)] converges in law when 
L — > oo to the t-independent Gaussian process (vq, w(r)) with the 2-point functions 

1 dk 



k d+2a > 

(l-e ifc - r )(l-e- ifc ' r ') dk 



(v v ) = D 2 

(Mr)MrV = D 2 J ^ (ar)J , 

(v w(r)) = 0. (7.2) 

Note the independence of v and w(r). It is then natural to conjecture that the following 
convergence in law takes place: 

L- a R(t) — ► v t, (7.3) 



L^oo 



L- 2 



R(t) - Jv(s,0)ds] — ► [w{v s)ds. (7.4) 

n J L^oo J 



In the frozen case or if a < (3 or a > (3 > a(l — a), the limits describe the leading 

terms in the single trajectory statistics for large L. For (3 < a(l — a), one should also 

t 

take into account the term coming from J[v(s,0) — v(0, 0)] ds which is of order L a_/3 . 

o 

The dominant term of order L a in R(t) describes the ballistic motion with the random 
velocity of the largest scale eddies that sweep the Lagrangian particle along. In Appendix 
3 we prove convergence (7.3) in the frozen one-dimensional model with a — |. 

How does the presence of large scale L in Eulerian velocities influence the Lagrangian 
particle separation? Let us try to understand this in the one-dimensional frozen model. 
We shall consider two particle trajectories x(t) and x(t) + p starting at time zero at zero 
and po > 0, respectively, and we shall try to estimate the behavior of their separation 
p(t). First notice that p(t) > 0, i.e. the order of the particles on the line will never change. 



2 In principle, we should add the noise to the trajectory but it does not play any role on the 
scales that will be discussed. 
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For large L, the dominant events are when the velocities of the particles and at the points 
between them are all of the order L a and of the same sign during the time interval (0, t). 
Let us suppose that they are positive, see Fig. 9 (the case of negative velocities can be 
treated in a symmetric way). The crucial fact resulting from the one-dimensional geometry 
is the magic identity 

yrJ x(t) Vr/ 

The left hand side is the time At that the first particle takes to reach the initial position 
of the second one. The best way to understand the above identity is by releasing the 
second particle after the delay At so that both particles move subsequently together. 
The delay changes nothing in the movement of the second particle since the velocity field 
is frozen. The delayed particle will then be at position x(t) at time t (together with 
the first particle) and at position x(t) + p(t) at time t + At. But the right hand side of 
Eq. (7.5) is the time that the second particle takes to move from x(t) to x(t) + p(t). 




FIG. 9. Positions of two right-moving particles in Id frozen velocity 



Hence the identity which may be also proven more formally by noticing that the time 
derivative of its right hand side vanishes. Writing for large L 

x(t) = L a v t + 0(L a2 ), v(x(t)) = L a v + L a2 w(v t) + 0{L ai ) , (7.6) 

see relations (7.3) and (7.4), and anticipating that p{t) = 0(1), Eq. (7.5) may be approx- 
imated as 
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from which we infer that 

p(t)-po = L a2 - a p v^w(v t) + 0(L a3 ~ a ). (7.8) 

The process w(x) is the fBm with the 2-point function (5.28) (up to normalization). The 
precise conjecture would then assert the convergence in law 

L«d-«) [p (t) - p ] _> p . (7.9) 

Note that the above calculations indicate that not only a single particle motion, but also 
the separation of trajectories in the Eulerian frozen one-dimensional velocity ensemble are 
dominated by the scale L velocities, i.e. by the large eddy sweeping. The effect on the 
pair dispersion is, however, inverse to that on the single particle motion. Whereas the 
latter one becomes very fast for large L, the trajectory separation becomes essentially 
frozen to the initial value in a localization-type effect. It would be interesting to know if 
such localizing tendency persists in the more general Eulerian Gaussian ensembles (3.1). 



That the sweeping modifies the pair separation statistics for finite L may be seen in the 
following way. There is a competition between two types of contributions to the dynamics 
of the pair dispersion p(t). The first comes from the configurations where the velocity 
differences at distances of order p(t) are much smaller than the velocity of each particle. 
The second one from the opposite regime. The two contributions may be separated if we 
fix the initial velocity v(0) of the first particle, with i>(0) < -D^Po corresponding to 

1 /2 

the first regime and v(0) > D 2 Po t° the second one. Denote by Q(t; po, pi\v(0)) the 
conditional PDF of the exit times for fixed v(0). In particular, Q(t; p , pi|0) is the quasi- 
Lagrangian PDF studied in the L — > oo limit in previous Section. As long as p , pi <C L 
and v (0) D^ 2 L a , the conditional PDF should be approximately L-independent and, 
consequently, it should satisfy the scaling identity 

pQ(pt;p>o,P>|p a -^(0)) ^ Q(t;p o , Pl \v(0)) (7.10) 
for a = -r—. We infer that 

Q(t; po, PiKO)) ~ pr 1 Q(pr 1 t;l, 7 |Po^(0)). (7.11) 

where, as usually, 7 = ^- Deep in the regime t>(0) < -D^Po the PDF Q(t; p , pi|i>(0)) 
is then essentially quasi-Lagrangian. As for the opposite regime, we may use the magic 
formula (7.5) with p(t) = 7p . Deep in that regime, the fluctuations of v(p) in both 
integrals are small and Eq. (7.5) reduces to the approximate identity 

v(x(t)) ~ 7«(0) (7.12) 
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from which p dropped out and which states that the exit time t is the first time when 
the velocity on the trajectory of the first particle reaches the value 7^(0). In the scaling 
regime, we obtain then 

Q(t;po, Pl \v(0)) - ^(0) 1 - 1/Q Q SC (^(0) 1 - 1 ^ ;7 ) (7.13) 
which is consistent with (7.11) in the crossover region v(0) = 0(p a ) ) Even for v(0) > 

1 /2 

D 2 L a where the scaling breaks, the po-independence of Q(t; p , pi\v(0)) persist so that 
the contribution of the region v(0) > D^ 2 p^ to the moments of exit time is approximately 
Po-independent for fixed 7. On the other hand, the contribution of the quasi-Lagrangian 
regime v(0) < D^ 2 pQ to the n th -moment is approximately proportional to 

J dv(0) jt n Q(t;p on p o \v(0))dt ~ Dl 12 p a Q +n ^- a) jt n Q(t;l n \0)dt. (7.14) 


It dominates for small po ii n < —j^- Altogether, we then expect that in the frozen 
one-dimensional Eulerian model and for small p , 

1 r, , v r x . , > i a + nil — a) for n < —-r— , „ „s 

(t n l {t<oo} ) = 0(jfr) with (n = 1 ' for n >Jj? (7-15) 

i.e. again a bifractal situation. The prediction seems to be confirmed, at least for large 
|n|, by numerical simulations, see Fig. 10. 




FIG. 10. Scaling exponents of exit-time moments in Eulerian Id simulations for two different 
box sizes. Illustration of the sweeping effects 
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Appendix A 



We shall establish here the estimate (5.9) on the average time, given by Eq. (5.8), that 
the trajectory takes to reach the first zero of w(p) between p and p 2 . To this end, let 
us note that for 8 > 0, 

OO 2 2 S 2 

f Q-^^w^w^ < ^ -3^ < ^ ( 6-4(^1 < y/2(AS) 8 B- 8 4g • (A.l) 



Employing this bound for < 5 < 1/2 and extending the integral over p in (5.8) to 
infinity with the use of the identity 

OO 2 2 

_ (mp —w) _ (w +w) 



J2ir(p-p ) 
po 



/ / e 2(p- P0 ) _ e 2(p- P0 ) g _ Wq _|_ w _ | Wo _ w \ ^ (A. 2) 

we obtain 



-5 „<5 f dw o [ (w +w-\w —w\)dw 







,,1+25 



v/tt <5(1— 2<5) 

The minimization over 5 gives the inequality (5.9). 



(2tf)«e-«r(l-fl 1/2-6 S ( A o\ 



Appendix B 

This appendix is devoted to the spectral analysis of the operator /C_ given by Eq. (5.21) 
with the attractive potential, pertaining to the long-time behavior of the d — 1, a — \ 
exit times. The two eigen-solutions of /C_ corresponding to an eigenvalue A may be 
expressed by the Whittaker functions 

ip x (w) = M i i(2 v /= 2Au'), Ww) = W i i(2 v /= 2Aw). (B.l) 

The spectrum of /C_ on the positive half-line and with the Dirichlet boundary condition 
at the origin is composed of the half-line [0, oo) (continuous spectrum) and of discrete 
negative eigenvalues E n = — for n — 1,2,... The eigenfunctions in the spectrum are 
ip E (w) where in (B.l) for E > we choose the square root with positive imaginary part 
and for E = E n the positive one. These functions vanish at zero. They are imaginary and 
oscillating at infinity for E > 0. For E = E n , they are real and decaying exponentially. 
In the latter case, the two eigen-solutions (B.l) become proportional and may be expressed 
by the Laguerre polynomials, similarly as for the three-dimensional Schrodinger operator 
in the attractive Coulomb potential: 

M nl C-w) = ^W n ,( 2 -w) 

2 v n ' n\ "■> 2 v n ' 
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= ^»e-i-l4. l( |«,) = Lyn*TL^ e -,)\ (B .2) 
The resolvent kernel of /C_ takes the form 

rr \Wf„, ,„ ^ 2 J^K)^W for w <wi, m v 

W Iv?a(w ) ^a(wi) for iy >iui, 

with the Wronskian 

W = <p x (w)d^ x (w)-^w)d w <p x (w) = 2 v /= 2A/ r ( 1 -7^)> ( B - 4 ) 

where in the expression for the resolvent the square roots are taken positive for A suffi- 
ciently negative and continued analytically to the other values of A outside the spectrum. 
The discrete eigenvalues E n appear as poles in the right hand side of (B.3) with the residue 

-f^ij Wo )tN (B.5) 

originating in the zeros of the Wronskian. Along the positive axis of A, the right hand 
side of (B.3) has a cut 



^(l -e ^e) 1 ^ e {w )^ e {w 1 ). (B.6) 
It follows that the spectral density of /C_ has the form 



00 2. v_1 



u(E) = ^2^-2E n 5(E-E n ) + ^(l-e-v*)" (B.7) 



n=l 

and that 



n=1 

00 00 

- / dwi / e-H 2 (pi-po)^ M i 1 (2i v ^E u , ) M 1 i(2iv / 2Ewi) d£. (B.8) 

Substituting this expression to Eq. (5.20), one can see that the contribution of the ground 
state of K- dominates for u = — i\u\ and large \u\ so that 

<e»'l (tol > = ^ e ^F^ (l + 0(\,n . (B.9) 



Appendix C 

As an illustration to Sec. 7, we shall prove here the convergence (7.3) for the one- 
dimensional frozen case of the Gaussian ensemble (3.1) of Eulerian velocities with a = \. 
Using the scaling properties of the frozen velocities, both x(t) and vq may be realized on 
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the same probability space corresponding to the velocity process v(x) with L — 1. This 
is done by setting 

x(t) = Lx(L~h), v = v(0), (C.l) 

where x(t) is the Lagrangian trajectory in the field v(x) such that x(0) = 0. We shall 
prove the (stronger) convergence (3.1) in the L 2 -norm on the probability space of v: 

([L*x(L-h) - v(0)t] 2 ) — > 0. (C.2) 

L— >oo 

If v(0) > then x(t) > and, symmetrically, if v(0) < then x(t) < 0. The contribu- 
tions of the two cases to the expectation (C.2) are equal so that we may restrict ourselves 
to the case v(0) > 0. It will be more convenient to estimate the expectations of the exit 
time t(x) of x through x > related to x(t) by the identity 

1 {x(t)>x} = 1 {t(x)<t} ■ ( C - 3 ) 

~ X H 

Since t(x) = J if v > on [0, x) and is infinite otherwise, we have easy bounds 

l{Vmi n >x/t} ^ l{t(x)<t} - ^{V^x/t} (C-4) 

X 

with v m { n being the minimum of v on the interval [0, x] and v av = - J v its average 

x o 

value. Now, with the use of the identity (C.3) and integration by parts, the L 2 -norm on 
the left hand side of (C.2) may be rewritten as 

oo oo 

V^V-Wo^ - 4t /<V-wi.> 5(0)><fa + t2<5(0)2) - (a5) 



From the explicit expressions for the Gaussian field expectations, it follows that, for t? min 
and -y av standing now for the minimum and the mean of v over the interval [0, L~^x], 

2™ (hv^x/t}) = J™ (hv^x/t}) = (l{*(0)>x/t})> ( C - 6 ) 

so that also 

&,<W*x)<L-i t }> = (lRo)>./ t} > (C7) 

and similarly with the insertion of £>(0). It is also easy to show a uniform in L bound 
(1|~ . >x / t y) < e~ c ( x /^ . From the Dominated Convergence Theorem, the limit of (C.5) 
is then equal to the expression 

oo oo 

^ f x (l {Zmx/t} ) dx - Atj(l { ~ mx/t} v(0))dx + t 2 (v(0) 2 } (C.8) 



which vanishes in a Gaussian ensemble. Generalization of this proof to the case with 
a 7^ does not pose much problem. 
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